Spatio-temporal Fermionization of Strongly Interacting ID Bosons 
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Building on the recent experimental achievements obtained with scanning electron microscopy 
on ultracold atoms, we study one-dimensional Bose gases in the crossover between the weakly 
(quasi-condensate) and the strongly interacting (Tonks-Girardeau) regime. We measure the tem- 
poral two-particle correlation function and compare it with calculations performed using the Time 
Evolving Block Decimation algorithm. More pronounced antibunching is observed when entering 
the more strongly interacting regime. Even though this mimics the onset of a fermionic behavior, 
we highlight that the exact and simple duality between ID bosons and fermions does not hold when 
such dynamical response is probed. The onset of fermionization is also reflected in the density dis- 
tribution, which we measure in situ to extract the relevant parameters and to identify the different 
regimes. Our results show agreement between experiment and theory and give new insight into the 
dynamics of strongly correlated many-body systems. 
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Interactions between particles are of particular im- 
portance in one-dimensional (ID) systems as they lead 
to strong quantum correlations. For ultracold bosonic 
atoms, reduced dimensionality and control of interactions 
can be achieved experimentally by quantum optics tools. 
This has lead to the observation of the strong-interaction 
regime of bosonic atoms []], One-dimensional ultra- 
cold gases offer moreover a unique test bench for theory, 
as they are among the few many-body systems which 
can be described on the basis of integrable hamiltonians 
such as the Lieb-Liniger (LL) model [j, Q. While the 
integrability provides exact benchmarks for a compari- 
son with experiments, important quantum characteristics 
such as higher-order correlation functions remain a chal- 
lenge 0, Their measurement would provide impor- 
tant information about the many-body system beyond 
that accessible from density profiles in coordinate or mo- 
mentum space [7]. Only local second and third order 
correlation functions, g( 2 )(0,0) and g( 3 )(0,0), have been 
experimentally investigated in ID via indirect diagnos- 
tics Of even larger interest are temporal corre- 
lations as they probe the nature of excitations, which 
goes beyond the characterization of quantum properties 
of the ground or thermal state of the system. In order 
to directly access correlation functions in situ, spatially 
resolved single atom sensitive detection methods are well 
suited [ll], [l2[. In this work, we use scanning electron 
microscopy to study one-dimensional tubes of ultracold 
bosonic atoms in the crossover between the weakly and 
the strongly interacting regime. We characterize the ID 
systems by performing in situ measurements of the spa- 
tial density distribution, which allows for the determina- 
tion of the relevant parameters and the identification of 
the different interaction regimes. The complete temporal 
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is then measured for £ = 0. Here \I/ are the bosonic field 
operators and (...) indicates the quantum mechanical av- 
erage. The results are finally analysed on the basis of the 
LL model, solved by numerical methods. 

The usual classification of the ID regimes for trapped 
Bose gases, under conditions that justify the local den- 
sity approximation (LDA), is based on the dimensionless 
interaction parameter j(x) at the trap center and on the 
temperature T [H, Q. The space dependent LL pa- 
rameter is defined as j(x) = mg/h 2 p(x), with m the 
mass of the particle, p(x) the density and g ~ 2ahu r the 
ID coupling constant for \a\ < a r , where a is the three- 
dimensional s-wave scattering length, uj r the frequency 
of the radial harmonic confinement and a r = ^/h/muj r . 
For 7(0) <C 1 and T below the degeneracy temperature 



Td = H 2 p 2 (0)/2mkB 14], with ks the Boltzmann con- 



stant, a weakly interacting quasi-condensate phase is pre- 
dicted. The spatial density profile in a harmonic trap is 
expected to be well described by a Thomas-Fermi (TF) 
parabola [13] and # (2) (0,0) & 1 [14]. In the opposite 
limit, 7(0) ^> 1, the ID gas is strongly interacting and 
approximates a Tonks-Girardeau (TG) gas. The density 
profile is a square root of a parabola and g( 2 )(0,0) <C 1, 
which indicates an effective reduction in the overlap of 
the particle wavefunctions, resembling the fermionic ex- 
clusion principle. 

In our experiment, a Bose-Einstein condensate (BEC) 
of about 8 x 10 4 87 Rb atoms is created in the optical 
dipole trap realized by a focused CO2 laser beam. Once 
the BEC is produced, its final atom number is accurately 
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FIG. 1. (Color online) Schematic of the experiment. The 
2D blue-detuned lattice realizes a system of parallel ID tubes 
of ultracold bosons. When the electrons (blue line) collide 
with the atoms they lead to the production of ions which are 
subsequently detected. The inset shows a zoom on the sum of 
100 images of the ID tubes. The integrated density along the 
x direction is also presented. In this picture the total atom 
number is ~ 10 4 , the current of the electron beam is 20 nA 
and its FWHM=120(10) nm. 

adjusted by scanning the electron beam on the outer part 
of the cloud. In this way we control the atom number 
and we selectively discard the warmer particles, effec- 
tively reducing the temperature of the sample (to less 
than 10 nK). For the experiments reported here, we pre- 
pared two different sets of BEC samples with 10(2) x 10 3 
(10k) and 60(5) x 10 3 (60k) atoms respectively. In or- 
der to create the ID atomic tubes, we adiabatically su- 
perimposed to the CO2 dipole trap a two-dimensional 
(2D) blue-detuned optical lattice, realized by a pair of 
retroreflected laser beams with wavelength 774 nm and 
waist 630 jira (see Fig. [1]). The final frequencies in the 
tubes for the two sets of measurements are co> a /27r = 8 
Hz, u r /27r = 56 kHz (10k atoms) and u a /27r = 11 Hz, 
u r /27r = 40 kHz (60k atoms), where uo a is due to the ax- 
ial confinement of the dipole trap. The number of tubes 
created varies from ~ 100 to 280, with a maximum occu- 
pation number in the centre of about 170 and 500 atoms 
respectively. Changing the density of the initial three- 
dimensional system and the radial confinement cj r , we 
are thus able to experimentally access two different ID 
regimes in the crossover between the strongly and the 
weakly interacting limit: 7(0) ^2 in the central tube for 
the set with 10k atoms and 7(0) ^ 0.5 for the set with 
60k atoms. 

In order to characterize in detail our experimental ID 
systems, we image the entire cloud by scanning it with an 
electron beam of 6 keV energy, 60 nA current and 240(10) 
nm FWHM [15j. At the end of the scan the number of 
atoms remaining in the trap is probed by time-of-flight 
absorption imaging and the BEC is checked to be pure. 
For each set of measurements ~ 300 images of the cloud 
are summed. Each scanning line of the image along the re- 



direction corresponds to the sum of the density profiles 
of the ID tubes displaced in a vertical row along the 
direction of the electron beam, see Fig. [TJ 

For the two sets of measurements, the density pro- 
files of the central row of the images, corresponding to 
the sum of respectively ~ 8 and 11 tubes along the z- 
direction, are shown in Fig. [2j The results are first com- 
pared with the exact Yang- Yang (YY) solution, where 
the only free parameters are the total atom number and 
the temperature. We numerically solve the YY equations 
for a uniform ID gas, for a given temperature, and we 
then apply the LDA to account for the axial confinement 
V(x) = muj^x 2 /2 [14J. The resulting density profile is 
then obtained as a sum over the different tubes, where 
the number of tubes and atoms per tube are determined 
on the basis of the TF distribution of the initial BEC. 
The final curves show excellent agreement with the ex- 
perimental data. Notably this allows for an accurate de- 
termination of 7(0) for the different tubes and of the tem- 
perature of the system. For the first set of measurements 
(10k atoms) we find that 7(0) ranges from 2.2(2) to 13(1) 
in the different tubes along the z-direction, resulting in 
the weighted average 7(0) = 2.7(3). The fitting param- 
eters are N = 9.3(0.7) x 10 3 atoms and T = 9(1) nK. 
For the second set (60k atoms): 7(0) = 0.58(5) - 2.3(2) 
and 7(0) = 0.76(7) with N = 52(5) x 10 3 and T = 22(2) 
nK. For further confirmation of the extracted tempera- 
tures, we additionally measure T in two ways: with a 
gaussian fit of the density in the outermost tubes along 
the ^/-direction and by calculating the average degener- 
acy temperature at which the density profile of the outer 
tubes starts to match with a normal distribution, when 
moving along the y-axis. We find: T = 10(2) nK and 
T = 8(1) nK for 10k atoms, T = 27(6) nK and T = 28(5) 
nK for 60k atoms. In Fig. [5(a)-(b), we also show numer- 
ical solutions for the LL theory at T = with LDA, 
calculated for the same 7(0) values extracted from the 
YY fit. The two models agree well, showing only a small 
discrepancy at the wings. This is due to the fact that the 
measured temperatures T are below the weighted aver- 
age of Trover the tubes, i.e. T [% = 0.7(1) (10k atoms) 
and T/% = 0.3(1) (60k atoms). To better clarify the 
regime of our measurements, in Fig. [2(c)- (d) we compare 
the profiles obtained with the LL model with those calcu- 
lated by solving the Gross-Pitaevskii equation in the TF 
approximation and by solving the TG model, with the 
same atom number [13[. For the lower 7, Fig. [2(d), the 
actual density profile is very close to the TF, while the 
TG distribution is significantly distant: the system lies 
in the weakly interacting TF limit. Increasing 7 above 
1, Fig. [2(c), the experimental data and the exact ID 
theory start to significantly deviate from the TF distri- 
bution. The TG profile, conversely, becomes closer: an 
intermediate regime towards the TG limit and the com- 
plete fermionization has been entered. 

When we proceed to the measurement of </ 2 )(t), the 
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FIG. 2. (Color online) In situ spatial density distribution summed over 8 and 11 tubes for the measurement sets with 10k 
atoms, (a) and (c), and 60k atoms, (b) and (d), respectively. Data (red dots) are the result of the sum of ~ 300 pictures. In 
(a) and (b) the green solid lines are fits using YY theory and the dashed black lines are T = LL profiles. In (c) and (d) 
the T = LL profiles are compared with those calculated with the TF (magenta dotted line) and with the TG theory (blue 
dashed-dotted line). 



differences between the two experimental situations out- 
lined above become even more evident. The observation 
of interaction induced antibunching in the two-particle 
correlation function is the most direct indication of a 
fermionic behaviour emerging in an ensemble of bosonic 
particles. The experimental procedure is similar to the 
one described in our earlier work Once the ID sys- 
tems are prepared, we investigate them by focussing the 
electron beam at the center of the cloud. What we mea- 
sure is the average correlation function of two ionization 

(2) 

events gkv (t) obtained by probing simultaneously 8 or 11 
different tubes. For these measurements we set the elec- 
tron beam current to 20 nA and FWHM= 120(10) nm. 
The spatial resolution, set by the FWHM of the electron 
beam, is smaller than the minimum estimated spatial cor- 
relation length for the probed systems, i.e. 440 nm (10k 

(2) 

atoms) and 170 nm (60k atoms). The gkv (t) is measured 
on the ion signal collected in 50 ms. This is the longest 
holding time at which the system is still a pure BEC. We 
compute the correlation functions over about 12000 (10k 
atoms) and 2000 (60k atoms) repetitions of the experi- 
ment. In Fig.[3jb) we report the details of the measured 

(2) 

gkv (t) for the set with 10k atoms, where antibunching is 

(2) 

visible with an amplitude 1 — gkv (t = 0) = 0.05(1). For 
r > the shape of the signal suggests the existence of 
two timescales on which the correlations relax: a short 
one (~ 100 /is) and a longer one (~ 500 /is). We first dis- 
cuss the longer timescale and, to this purpose, we show 
the data with 80 /is bin time in Fig.[3fc)-(d). From these 
data, using an exponential fitting, we extract the corre- 
lation times t c = 240(80) /is (10k) and 100(40) /is (60k), 
defined as the times corresponding to the 1/e reduction 
of l-g£\r). 

The g^ 2 \^r) function in the ground state can be the- 
oretically calculated for each single tube, by simulating 



the dynamics of the system using exact numerical meth- 
ods. We note that the use of a T — theory is justified 
in this case since the majority of the investigated tubes 
are below [16]. In order to perform the simulations, 
we discretize the LL model [17[. This lattice problem 
can then be treated using the Time Evolving Block Dec- 
imation (TEBD) algorithm [l8|, which is one of the ex- 
tensions of the Density Matrix Renormalization Group 
to time evolution [19j an d is based on a Matrix Product 
State (MPS) ansatz [20[. We perform the calculation on 
a gas of TV = 25 particles confined in a harmonic poten- 
tial, g^i^r) is evaluated in the center of the tube at 
xo = 0, according to the experiment. The numerical pro- 
cedure consists of first calculating the interacting ground 
state |0) at fixed particle number N • We then apply 
the annihilator at xq = 0, which is straightforward in 
the discretized MPS representation, resulting in a N — 1 
particle state |^o) = ^|0). This state is finally propa- 
gated in time. The spatio-temporal correlation function 
<7^(£> r ) = P~ 2 (^r\p( x )\^r) is then the expectation 
value of the local density operator p(x) = &(x)^(x) in 
the time-evolved state |\I/ r ) = exp(—iHr/h)\^fo), where 
H is the Hamiltonian of the system. The finite size of 
the system turns out to have no influence because all the 
excitations, travelling with the speed of sound, reach the 
edges of the cloud only at a time when the antibunch- 
ing structure in g^ 2 \r) has already decayed. In addi- 
tion we take into account the finite width of the elec- 
tron beam by averaging over its Gaussian profile W(x'): 

9 {2 \ r ) = J-oc dx ' W ( x ') (2) (*'>t)> see Fig. Ha). As 
the experiment measures several independent tubes at 

(2) 

the same time, we calculate gkv (t) by averaging over the 
g( 2 \r) corresponding to each tube [2l|. The uncorre- 
cted signal coming from different independent tubes is 
responsible for the reduction of the antibunching ampli- 
tude from the value of 1 - <? (2) (0) ~ 0.60 (0.40) for the 
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FIG. 3. (Color online) (a): g (2) (r) calculated using the TEBD 
algorithm together with the respective Fermionic counterpart 
for the single central tube, (b): Measured gj$ (r) for 10k 
atoms. The solid line is obtained from TEBD numerical cal- 
culations of the independent tubes [2lj] . The theoretical curve 
has no free parameters, since they are determined from the 
analysis of the density profiles, (c)-(d): Measured g$ (t) 
with binning time of 80 jjls for 10k and 60k. The solid lines 
are exponential fits to the data, (e): Correlation time t c as a 
function of 7. The dotted line are bare TEBD calculations, 
the solid line includes averaging over the electron beam size 
and the dashed line shows the fermionic theory, (f): Ampli- 
tude as a function of 7. The uniformly filled and patterned 
areas are the theoretical calculations according to 21]. The 
extension of the area is due to an uncertainty ±1 on the num- 
ber of tubes. 



single tube to the global calculated 1 — g& v (0) — 0.07 
(0.03) in the systems with 10k (60k) atoms. Concerning 
the longer timescale, the agreement of the averaged the- 
oretical curves with the experimental data in Fig. [3fb) 
is fairly good. For a quantitative comparison, we show 
in Fig. [3^e) and (f) the experimentally derived correla- 
tion times and amplitudes together with their theoretical 
counterparts. The correlation time t c rescaled by 1/p 2 is 
independent on p and is strongly depending on 7 in the 
regime of our measurements: increasing the interactions 
above 7 = 1 we observe a significant decrease of t c . This 

behaviour is in good agreement with our experimental 

(2) 

findings. The interpretation of 1 — #av (0) as a function 
of 7 in the regime that we consider is mainly complicated 
by the average over a different number of tubes. This 
effect is indeed contributing to the increase of the ampli- 
tude for the higher value of 7. Anyway the comparison 
between the experimental and theoretical amplitudes is 



compatible within the error bars. We note that finite 
temperature effects cannot have a significant role in the 
reduction of the amplitud e |16j |. 

It is well known [22l424j that ID bosons can be 
mapped to interacting fermions for arbitrary interaction 
strengths. As a consequence the density distribution 
and density-density correlations of bosons are identical to 
that of the corresponding fermions. Conversely, this sim- 
ple Bose-Fermi duality breaks down when the temporal 
pair correlation function is probed: the theoretical curves 
in Fig. [3]^a) and (e) show indeed a notable difference be- 
tween the fermionic and the bosonic case, which is con- 
firmed by the experimental data. The observable mea- 
sured and calculated in this work, i.e. r), is fun- 
damentally different from the dynamical density-density 
correlations previously derived in Ref. [25| because the 
bosonic field operators do not commute at different times. 
This is also the reason for Bose-Fermi duality to fail. We 
finally discuss the shorter timescale visible in Fig. [3f b), 
which corresponds to an additional modulation of the 
signal with frequency ~ 5 kHz. The energy scale intro- 
duced by the electron beam can be roughly estimated 
to be ft 2 /2raFWHM 2 -/ix4 kHz, suggesting it to be 
the origin of such an oscillation. We note that a similar 
modulation is not visible, within the error bars, in the 
weakly interacting case and was never observed in a 3- 
dimensional degenerate sample This suggests this 
behaviour could be associated to the presence of strong 
correlations between the particles when 7 > 1. 

In summary, we have performed in situ high-resolution 
measurements of the density profiles of few ID tubes in 
the crossover between the weakly (7 ±2 0.76) and the 
more strongly (7 — 2.7) interacting regimes. We find ex- 
cellent agreement in the comparison with a model based 
on the YY solutions and LDA, allowing us to extract the 
7(0) values of the tubes, the global atom number and the 
temperature. Comparison with the T = LL, TF and 
TG models allows us to clearly define the regimes of the 
measurements and to circumscribe the role of tempera- 
ture. We measured the temporal two-particle correlation 
function, observing genuine antibunching in correspon- 
dence to the entering of the more strongly interacting 
regime, characterized by evident fermionization of the 
bosonic particles already at intermediate values of 7. We 
have numerically calculated <7^(£,r) and we find fairly 
good agreement with the experimental data, showing the 
exact duality following the Bose-Fermi mapping does not 
hold for this observable. This work paves the way to 
the study of complex dynamics as those resulting from 
quenches j26[ and to the investigation of thermalization 
processes in ID [I?}. 
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